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1. Introduction 


Although the literature on nonlinear magnetic effects is extensive, there is little quantitative 
information available on multi-frequency excitation of ferromagnetic systems, even at the level 
of perturbation theory. The calculations presented in this technical report constitute an attempt to 
remedy this situation. Because it is in the nature of remote sensing to involve small probe fields, 
the analysis needed is an ideal use of perturbation theory, cumbersome though it may be. It is 
hoped that their derivation will help to lay the groundwork for evaluating the usefulness of 
ferromagnetic resonance (FMR)-based detection of the nonlinear excitation of magnetic objects 
in a true remote-sensing environment, 1.e., where large magnetic fields of the sort used in 
laboratories or medical equipment are unavailable. 


2. The Landau-Lifshits and Landau-Lifshits-Gilbert Equations 


In order to estimate the power radiated by FMR excitation of a finite magnetized body, it is 
necessary first to describe the dynamics of the magnetic dipole moment per unit volume M (i.e., 
the magnetization) within the body in the presence of an external time-dependent magnetic field. 
This dynamic problem is governed by the Landau-Lifshits equation (/): 


CUE eerie — ——__ ff x —— 
Ot eff IM | Ct 
St (1) 


where M st 1S the static magnetization of the material, y 1s the gyromagnetic ratio, o is the 


damping constant, andH _,, is an effective magnetic field consisting of the DC anisotropy field 


eff 
H_, that fixes the direction of M st 1n the material, the depolarization field due to the body’s 


shape, and the external time-dependent applied field. Note that in the absence of an AC magnetic 
field there is no depolarization field, and so M ap H q = 0 under DC conditions, 1.e., there is no 


DC torque. The first term is the torque exerted by the effective magnetic field H__,, while the 


eff 

second term gives rise to damping via eddy currents. The form of this equation implies that 

- 0M <2. — 

M: a =0> M| is constant in time. 
if 

Let us solve this problem approximately using perturbation theory. It is more convenient to first 


put this equation in Landau-Lifshits-Gilbert (LLG) form (2): 1f we define the torque vector 


= _ = : 2 a 

Q=yMxH eff and a dimensionless vector damping rate Q = ——— M , equation | can be 
Mr 

rewritten as follows: 


M es. oe et 
oa =(1+07) (Q-Ox4) (2) 
Ot 

Assume a coordinate system with the z-axis along the wire axis. We start by defining a small 
perturbing radio frequency (RF) magnetic field h, which generates an RF magnetization m that 
is small compared to M st: Lhe decay dynamics of this magnetization is strongly affected by the 
conservation of length of the total magnetization vector, whose endpoint is constrained to lie on a 


spherical surface at all times. Let us further specify the perturbed magnetization and effective 
magnetic field in the LLG equation as follows: 


M=M,,+m (3) 
and 
H og = Ha th—Nin, (4) 
where My xH, =0 and 
N 0 O 
N=|0 N 0 (5) 
0 0 N, 


is the depolarization tensor. In this report, we consider a long wire parallel to the z-axis, for 
which N =1/2 and N z © Q. Following Antonenko et al. (3), we define a static susceptibility 
Ms 


st = = . Then in rectangular coordinates with M st and H q along the positive z — axis, Le., 
a 


My =|M [2 and H, = Hgz, we obtain 


a =h,—-Nm, 


: (6) 


Hy =A6¢ |Ms|2> ee =hy —~Nm 
ie =h; +X5¢ [Mel 


A A A 


x y Z 


=> O=7(M yy +1h)x Hoy = my, my IM ,,|+mz 
re 4 lie 
h, -Nm, hy—Nmy he+X 4 M op 
Ae, ee Py (7) 
x -|M ,,|hy + Nis +N M ,,|my +myh,—-mzhy+Nm,my 
allay 1 \ag ~ 
= +5||M 5e|hx —(xs¢ + NM 5,]my +mzhy—myh, — Nmzm, | 
+2|myhy -myh, | 
Introduce a scaled dimensionless magnetization g¢ such that m= M st E . Since M st = M rr | Za 
the length constraint on M becomes 
S92 2s i ~ 42 wo! aan g) 9) 9) 
M =|M,,+m =|M,,| =>|2+g/ =(l+g7) +2x +gy° =1 a 
2 2 
=> 82 =—I+\l-8x — By 
Let Q=a(2+8). Then 
M (1402) (B-Gxai (9) 
St Ot 
and 
a hy +(X5 SEI, 
Q=y|M |} + 5 hy -( as 1 +N) M o|8x+8chy- ghz -N |W gece (10) 


+2| gxhy —gyhy | 


Defining the constants A = (x +N A 


O= IM é | N lets us write the torque vector as 


&| hy +Agy +gyhz—g7hy+Ogzgy | 
Q =7|M s,| +$[hy-Agy, +gzhy—gyhz -Ogzg, | (11) 
+2 [gyhy -gyhy | 


The right side of the dynamic equation 2 can be treated as a “source”’ S: 


Sx 


x * z Sy =(Q-QxQ) =O, -a(gyQz—-[1+27]Qy)=|M, 
QxQ=algy gy 1+g,;>4Sy =(Q-QxQ) =Qy -a([1+ g7]Q,-8,Q,)=|M,|Sy (12) 
ety Oy ne Sz =(Q-QxQ)_ =Q, -0(g,Qy -gyQy)=|M o,|8- 


where 
“hy +Agy +gyhz —Szhy + O8 78 y 


Sx =¥4 78 y (exh y -gyh,} 

+a[1+g-](hy —Agy +gzhy—gyhz -@g 72) 
hy —Agy +8zhy —gyhz —Og 78x 
-a[l+g](-hy +Agy+gyhz—g-hy +Og-2y) 


N> 
[| 
<2 


TOE x (exhy -gyh,} 
gyhy —gyhy 
S, =7)—agx (hy, -Agy, +g7h;,—g,h, -Og-g,) 


+ogy(-hy +tAgy+gyh;z—gzhy +Og-gy) 


Note that the vector S$ has the dimensions of frequency. 


It is advantageous to write the transverse system variables in vector form: 


Then, the equation for the transverse components of the source vector can be written as 


= —l., 
 =(1+02) S 
Ot 


: 2 
[ley all+g-) —[l+gz] er ae Bind ec Resta 
[i+g-| afit+e,|° A+h,+@g,-aG@H_ al[l+g,|[A+h,+0g, | 


a| 1-6? | Vi? ai-G?) A-O +h, +0 1-6? | -|A-O+h, +6 I-G? 0G eH | 


= H-y 


1-G? a 1-2 | A-@+h,+@yVI-G* -aG@H aVi-G?| A-O+h, +0 I-67 | 


(13) 


(14) 


(15) 


@ Uo) 


Here g,° +By* G? and G@H[) 2-GxH =xhy —Zyhy. There is also an equation for the 


longitudinal part (z-component): 


(1+07)<s, +ya(G-H-0G7)g, =y{G@H-oG-H+a[A+h,]G*| (17) 
[ 


where G: HO g xh, + gyhy , but the algebraic constraint on the length of the vector g¢ makes it 


redundant. 


Up to now, all the dynamic variables g,, g y> &z andh,, h y? h, have been real. In the spirit 


of the Holstein-Primakoff transformation (4), let us define the following (complex) “circular” 
basis for the transverse variables using the matrix 


— ., {8x +lg 
me b : 8x —'8y G 
— 
l 


C= . (18) 
es ~. (hy tihy H- 
H=CH= = ¢ 
hy ih y H-. 
In this notation, 
1 Xk 2 1 ok 2 
contests -tyheo(Geg" ete) 16-6 tome 
(19) 


= —(G He -GHe | a Im(G He | 


zi Z e Aa —B 
The matrices Z and R in equation 3 are both of the form T [ A, B| = a , which can be 
OL 
diagonalized by the similarity transformation defined by C: 
aoe ~ y {1 1 \(Aa —-B\1/fi 1 aA +iB 0 
CT|A,B|C * = | — = | (20) 
1 -i)\ B Aa)2i\1 -1 0 aA —1B 
leading to the following expression for the source vector in the circular basis: 
S S, +iS fae utes 
}- oS |=yEH - IG (21) 
S_ Sx — Sy 


where the matrices 


(22) 


and 
aJl-|gP A-06n, vot-(gt | 
0 
Aon, +0/1-|Gl -ain( 9H) 
Il = (23) 


afl (GP | A-o+n, voyi-igP 
+] Aon, +041-(GP -a1m(G"He) 


are both diagonal. Since the nonzero elements of these matrices are complex conjugates of one 
another, the system reduces to a single complex-valued equation 


1] 044, +O1-IGf ~aim(G"H1.) 


| 2 
where G and G * are “magnon” variables (5). Note that the quantities 4/1— IG | and 


Im(G" He are real. 


= rl-IgP ‘ay t-lgp + j H,(24) 


(I+07)“G+y 


A straightforward power-series expansion of this equation up to third order in G leads to the 


following perturbation equations: 
(a. JEG, +yAG, = 7H 
(a-i)<G, + AG, =-yhzG, a 
(1-1) 2G, +1aG, = (a +i)heGy 304+) Gi] He razG Me 
t 


Q+1 1 


H(ar+farJo)lgl G, 


The corresponding z-component equations are 
Z xk Lee 1 * * 
gz =-1+y1-|G| aly laG sleere G-5(G G,+G, G+ 


871 7 (26) 
asi I . 
> 1879 = 5% G 


573 ~ -=(G/'G2 as Gy Gi) 


Note that according to these expressions an incident field with only h, nonzero, 1.e., polarized 


along the direction of static magnetization, cannot generate a /inear response. This is in keeping 
with elementary considerations. 


3. First-order Solution 


To solve the equation 
\ Oo 
(a-) Gy + AG, = 7H; (27) 


we make the following decomposition: let 


Then 


9) 8x ~| 8x ~ {Sx DN ins ~ ~ a 
(I+-0 ale Joga B= ' J+(0 ]0,8 +0 p Mg = MS 


aero ba 
e-(| eee bl Sy §c— 2i\ 1. HL) Bes 

1 -i Sy Sc Sy _@l Sct |_ L(t ft \f Sey (28) 
Sy} Sat) Or le SIS 


_ 7 eT A 
= 1+a7)o,€ {tJ MC Met) "| 
§c— §c— Sc- 
ee eT a 
=(1+a7}o (f+ J EMC i fer) em i a 
t B 5 
Ec- Ec C— 


Sete, (Le. ele =e ige 3 
CMC = 
1 -z/\l a | —l-ia i 1 -l 
tol ia—l—-l+ia ia-—1l+l-ia ail 2ia — 2 ati @ 
«il iat+1—-1l-ia iat+1+1+ia) 2i ia +2 0 a-i 


rontehetetepteie) 


The following is complex notation for non-Hermitian quantities: 
S S S S : ; Soc. ae a 
rod : so orf 5 fin ia se jefe veton)( 5 fo ior] 
y ye ys Soe ys } 2 
Zan xs ior | Sxe ~ xs plot Sct tot | Sx plot 
9) Syco FiS ys Syco iS ys Sy4 Sy 


a ee ior , Sy— +S y— jor _ Sxc tiSxs +S yo —Sys _iot , Sxc —Sxs tS yo + Sys Jiot 
Sc- ys 


S48 Se +iSx5 —iS ye + Sys Sx¢ —iSx3 —iS ye —S 
Sxo—Sys +i] Sys +5 


(30) 


The following is the transient solution: 


—p f _ ‘ 

Op 8ce4(t)= Aye T + Bie Cag 
—p t _s : 

GAs (t) = A_e ey a aes ol 


= —1ot iot _ iY —10ot * it 
(2, +P, )&e+ =(P, io) Bye +(p, +i0)C,e =—1_(S,¢ +S_ e 


ioe 5 ee, 
4 a-i p, —iw ® p —i@(a—i) Op —O- 100 
—> 
= ay s 


a-ip, +i® Op +i@(a—i) Op t+ @Otiaw 


(2, + p_)&e- = (p_ — io B_e Mt + (p_ + io)C_e'™ See (Sze + Seo) 


ati 
eel p=)... yS_ 7 yS_ 
—  atip —i@ Op —i@(a +i) Op + O-iaw 
—> 
Bt se Sa" Sy" (31) 


a+ip_+i® Op +i@(a +i) Op —O+ 10 
For the boundary condition at t = 0, 


Bet (0) = 849 = Ay + Be + Cy > AL = 8 e439 ~ Be — Cy 


—pt KY - —pt co . —pt 
8ce+(t)=8 0° - ee be ae Bre 2g or aaa BEE eg ok 
Cr Op —@—ia@ Op +@+iaw 


Sy ] 
5 L8c+ Sc | 
—pt Se —pt S ae: 
GSR See, —-,_1 iot_, “- |,_1 + lot _, 
+ . c—0 
a+ip —io a+ip +i@ 
se *k *k *K *k *k 
pit S_ _ pit S —p Tt 
=e ae + L) : lot a + = gimt _, “+ 
Om O-—I Pp —1i10M Q—l D +10 


* * k kK * * 
—p tT S= —pit S = —pit 
Q—1 Py —10 Q—I zn +10 


S =f S_ 
go (1) > YOu 2 iot Y | plot 
Op, —O—i0M Or, +0+iA0 
B B 
Large t > 5 
S = S 
ee (1) > v = iot YOu | iwt 
Op +O-iA® Op —O+10® 
S cc oe 
Re YOu | SIOP «. Ds | plot 
pe Op — O10 Op + O+ 10M 
=> : = : 
y Ss _ ; 
fea YOu 7 jot ys | plot 
Op — 0-100 Op +O+iaw 


The linear source terms are 


10 


Q 


°x = nx = ox cos wt + Dxs neese Dxe + xs e tot Dxe ~ Hxs elt 
Sy hy hye hys 9) hye + ihyg hye — ihyg 
* * 
(ln n hy plot ean ae a in hy plot 
xk ~ ; *K 
hy hy Soe 1 -il hy hy 


* OK x (36) 
= nerty —i@t ey iot ,{ 44 )\ -ior | 4- | ior 
= le ite, <a Ze e + x le 

hy, —thy hy —thy A_ Ay 
—i0 
eG 
2 2 2 2 x x 
Hy = lyc” thyg” , Hy = hye thys = -i0, 
hy = Hye 
The two-frequency solution is 
8x t1Zy Gi 
= * 
8x —!8y G, 
Y —1Mt * 1Mct 
0, + [H S +H_“e % | 
fa i ai"? 


* 


= Wise + isy) otal 1 (hsx —iltsy ) {itet 


G1 Op — Wg — 10LWg Op + Wg + 100g 
= *K 
cos tiny). mop, 1x =ipy) fog 
Ip Op — Op —1AMp Op + Op +1AMp (37) 


4. Second-order Solution 


The solutions to 
0 


require the first-order solutions in the forms 


11 


10) ; 
Gia b * abt" abe? af ah ype - 


m 10) pf Z a HO pi oe 
Gia.b 7 ey ee cy pe De ° 
which gives the following expressions for the source terms: 
7 7 OT Opt { ! 
yhzG, = Re{ hg )s Reh. Gia + Gis 
iW t *K * +10 gt 
ae . Aad wates Eas ieee 
= ae tity se se e +h, e ae , : tet 
tT, Aye iy Ay e€ 
nv oH 62a sh rH, +h 1, Ho sum yy or He ait 
az at at az @ ~ a— az b+ b+° az b-  b- 
+10 ,.-t iw, t 
lO sym! . . dif b a‘ ‘s 
_y th I, i, ,¢ +h, Va Hy e +h, 1, Ap? led ps Ay 
2 & is * i eg i +10 di t 
a if He TO sym! 
re r gtd os a Pils Dy Hg e acre D4, e ure I, Ap é (40) 
x LOGE x x * LIM t x x x x +210,¢ 
sum b 
+h, | l iH, ,e +h A de +h, Day th, l, Ap e 
[1 row eo +(h ee) ot iOal 
az a+ at az 
h Vo Ho +h. TV, H, +h. VT. H . +h T..H 
a a Te b- = az at A ps b+" b+ 
eh 1M sympt ( kk K ko ‘ +iMsymt | (41) 
; +(4 Ty, Hp, +h, UH, Je +(h, Ty, Hy +h, Ta Ha- |e 
‘ é . —10 ae Z : o +1@ dif! 
(Ty Ay +h, Ls iA, Je +( hy Pa Hg tis Py Apa é 
x ok x +21@,t —2i@,t 
pss ~~ as 
where ® un Og tp and @ dif =@  —@) are the sum and difference frequencies. The basic 


equation for the second-order perturbation G 5 1S 
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( h x H eee +| ete) ot 2i@al 
h Vg Hg +h V, H, +h VT.H +h, T,,H 
TMqz a— Fa- T"p' p— bp az at at "bz ” b+ bt 
. O a —1O symt a ex *) +10 cmt 
(2-1) Go + ORG ar +(h Ty. My +h TH le +(h,, Ty, Ap +h, Va- Hae Je 
“3 . sy Oe i 3 7 HO) ee 
(a, Ty Ay tM rH le +( Pa A gq tel Ty Apa le 
Ok x +21@,t —2iM, t 
b b 
bs i ee LS ec eh 
Pu pie ee. oe ae, a ee 


Db OOF 


Pig he a, ee a Ee 
+r h Da Bee pat ge Git gal” _ 


Pat at azat 


bz b+ ob 
21M pt sk a a +21@ 1 
> Gy =-| 4 5, Vy Api? se Pe 
* * * 10 gyjyyt * * * * * * +10 gyyyt 
Ls aig a: (1 Ty Ay +h, r iA, a HD sym (1, DA li Vqg- Hq Je (42) 
. i" Be Oe é se re Bir 
ra Freee Hy eh Td. Je + Paces Hg +h Ty Apa le 
where 
1 1 
ae a 
2 Op = 2 Op t 20g — ia(2@zq ) 
1 1 
i ——— ro, =t—_________ a3) 
~ 2 Op F 20, ~ia(20, | ae Op + Osym — 12(Osym ) 
geen 
(ais ) 7 


Op Og? ~ io{ gp 


are new resonances associated with the various mixing products. For the z-component, we have 
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aoa al 6] 


9) we ; 2 
—iw_t sk * +1°0,f 100, t se ste TU, | 
Teydald aan <r SM es es Ppa Hy aa Hy e 
2 -i@pt * +i@t IO) t *  « T1Oyt 
+2 Re Tas Hase Fl = Aga € Ape ao Ay e 
OT sk sk +1O yt 2 IO, se se +101 2 
; Leda ce lj Hae + TAy aa ae Ay e 
a) = 
—iw_,t sk * +LI@_,ft LO t ste ste TIO t 
aR [Fy Mg cle cl te «re J Tot tle Fp e 
Z x (2 2 2 —210 jt 
gpa: eral | pea ¢ +2Re/Tyg Hg Tg4Hy4€ 
| HU, Ay, Pp Hy +2RelT, “H, T, Hye 
2 1 symt . : 10 ig 
12Re(Pg Ha Ty, Hpyé )+ Rel Hat, Hy e 
+10 ,.-t ; 
42 Rela pH De Oe ‘ 2Re(F A: rl Os age 
x x sf * +I t —10 t 
Re(T a, A a4 Dy _ Ay sum J=Re(Pa Hast, Hy_¢ sum 
+10 ,.-l 1M ,.-f 
dif *k *k dif 
Re Fa Aas re ea J=Rel Pastas, Ay. e 


2 x x |2 * *k —210,t 
Te aed pall el ges ge +2Re( Tg i; Pe Ro EE: 


a : 2RelT, H, T, H aa 
+ b-— ~ b- Do rps pc Teen pee 


+2Re([Tq_HaT H, +0 q,H,.0 


10 gy ppt 
b+ b+ _Hy_|e | 


b (44) 


x x x x aoa! 
+2 Re pda hs, Ay +l qyHasl,, Ay, le 


Two frequencies of special interest to remote sensing are the sum and difference frequencies. 
The magnetization amplitudes at these frequencies are given by 
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: * ; 
ms —1O sql ( H1O oy yyt 
Go com = | aan (4,5; Ay tl Ay Je +T sum—(#, 0), 4, +4, 0aHa-} ¢ 


8 2sum,x Re Go sum (45) 


52 sum, y IMG 5 sum 
ADeinit 
_ -Re([Ty_H ee DEE Pee ee 


29 sum,z a~" b+ b+ 
and 

ie = sty et ss a HO ae! 
Goaif | gpalbe ps A ale Je era hy, Va— Ha +h, ry Ap le 


Bodif x = RG aaig (46) 
Soait.y = ™G rai 
* 2 2 * tO aie 
8 a dif ,z =—Re gst Si Ay +P Aail yp, Ay le 
5. Third-order Solution 
The time-domain equation of the third-order system has a complicated source term: 
1 AIA If Me er 
6 ; -5 (30+ )/G,| Ht +a-Gy Ht 
(a -i)—G3 +OpG, =—— 
Ot aO+i} | 2 
+S (ad +[a+i]0)|G, G, -(a+i)hzG, (47) 


y 20 oO o 
—<—| 1+ S, + S, + A+@ |S, -—25 
L a Ll ati 2 (= 3 4 


To evaluate G3; we must first express the various products in multi-frequency form. This leads 


to the following linear inhomogeneous equation: 


15 
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+53 sa4b4¢ 
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+53 94 -b 4° 


which is easily solved: 
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05 °3304° 
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FE ahi 3. dnb” 
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—310 yt 
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\t . es (204 +0 

5 a \t P hee (20, — Og )t 
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43,548 tA pape 
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Note that the third-order source terms give rise to three distinct effects: frequency tripling, 
amplitude-dependence of the response functions at the drive frequencies, and intermodulation 
products. 
6. Grouping Third-order Magnetization Terms 
The only terms of interest here are the intermod terms, which are color-coded to make the 
algebra easier to follow. 
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7. Conclusion 


This report was written with the intention of archiving the results of lengthy calculations needed 
in order to evaluate the feasibility of remote detection of magnetic systems by FMR. More 
specifically, these calculations make it possible to calculate the induced magnetization within a 
ferromagnetic body illuminated by electromagnetic waves at levels of excitation high enough to 
give rise to a nonlinear response in the body and induce it to radiate at new frequencies that are 
detected at a remote receiving antenna. The results tabulated here contribute to the study of 
nonlinear radar and target signatures in general. 
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